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The elHptical instability of a rotating stratified fluid is examined in the regime of small Rossby 
number and order-one Burger number corresponding to rapid rotation and strong stratification. 
The Floquet problem describing the linear growth of disturbances to an unbounded, uniform- 
vorticity elliptical flow is solved using exponential asymptotics. The results demonstrate that 
the flow is unstable for arbitrarily strong rotation and stratification; in particular, both cyclonic 
and anticyclonic fiows are unstable. The instability is weak, however, with growth rates that 
are exponentially small in the Rossby number. The analytic expression obtained for the growth 
rate elucidates its dependence on the Burger number and on the eccentricity of the elliptical 
flow. It explains in particular the weakness of the instability of cyclonic flows, with growth 
rates that are only a small fraction of those obtained for the corresponding anticyclonic flows. 
The asymptotic results are confirmed by numerical solutions of Floquet problem. 

1 Introduction 

Elliptical instability, the three-dimensional instability of two-dimensional fiows with elliptical 
streamlines, has been the focus of a great deal of research activity. The review by Kerswell— 
discusses the main results up to 2002 and emphasises the relevance of elliptical instability to a 
broad range of applications. One of these is the instability of two-dimensional vortices that are 
deformed elliptically by a large-scale strain flow. This is especially important for the dynamics 
of the atmosphere and ocean, since this is characterised by an abundance of vortices that are 
deformed through either mutual interactions or the effect of large-scale flows. In this context, 
however, the planetary rotation and density stratification need to be taken into account. 

Rotation and stratification clearly exert a strong influence on elliptical instability: since 
this stems from the parametric resonance between the periodic fluctuations associated with 
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the elliptical motion and the free waves supported by the flow, the dispersion relation of these 
waves is critical. In the presence of rotation and stratification, the waves are inertia-gravity 
waves whose frequency is bounded from below by the minimum of the Coriolis parameter / and 
Brunt- Vaisala frequency N. As a consequence, a vortex of fixed vorticity ceases to be unstable 
by the subharmonic instability responsible for the simplest form of elliptical instability when / 
and exceeds a certain threshold. As these parameters increase further, instabilities are limited 
to resonances of higher and higher order, leading to decreasing growth rates. This was clearly 
demonstrated by Miyazaki^ on the basis of numerical solutions of the Floquet problem that 
models elliptic instability. Further numerical results were obtained by Mc Williams and Yavneh— 
who concentrated on the regime of rapid rotation and strong stratification with N > f most 
relevant to the atmosphere and ocean. Their broad motivation was the role that instabilities play 
in the generation of inertia-gravity-wave-like motion, and the resulting breakdown of the nearly 
geostrophic and hydrostatic balance that is typical of much of the atmosphere and ocean. The 
present paper shares the same motivation. It re-examines the elliptical instability of a rotating 
stratified fiuid and derives explicit analytical results in the limit of fast rotation f ^ Q and 
strong stratification N ^ il., where denotes the (relative) vorticity of the fiow. 

Several recent papeis^^'^ have demonstrated in specific examples that instabilities of well- 
balanced basic fiows to inertia-gravity-wave perturbations (or perturbations related to similar 
fast waves) have growth rates that are exponentially small in the Rossby number, here pro- 
portional to <^ ll^ Theoretical arguments^i^ indicate that this is a generic property, and 
the elliptical instability examined in this paper is no exception. In this case, the exponential 
smallness can be roughly understood by noting that, in the manner typical of parametric insta- 
bilities^, the growth rates of the elliptical instability can be expected to be proportional to 0", 
where n is the order of the resonance. Since, as pointed out by Miyazaki^ , the minimum n is of 
the order f /il. (for N > f), this leads to the conclusion that growth rates are beyond all orders 
in the Rossby number. To go further than this rough argument and provide an estimate for the 
growth rate requires the exponential-asymptotics analysis of the Floquet problem relevant to 
the elliptical instability. We carry out this analysis and, rather than relying on general asymp- 
totic results for Hill's equations directly relate the growth of solutions to the occurrence of a 
Stokes phenomenon— which we capture using a combination of WKB expansion and matched 
asymptotics in complex time. The analytical results are confirmed by the numerical solutions 
of the Floquet problem. 



^Note that we use the convential form of Rossby number rather than its inverse as used in Ref.- 
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One of the issues which our treatment clarifies is the difi'erence between the instabihty 
of cyclonic and anticyclonic vortices. Cyclones have been recognised as less unstable than 
anticyclones, to the extent that McWilliams and Yavneh'^ considered only the instability of the 
latter. We show that the cyclones are in fact linearly unstable, with growth rates that have 
the same exponential dependence on the Rossby number as the corresponding anticyclones but 
differ by a factor which, although formally of order one, turns out to be numerically very small. 

The plan of this paper is as follows. In section 2, we formulate the problem of elliptical 
instability in a rotating stratified fluid modelled using the Boussinesq approximation. We use 
the simplest instance of elliptical instability, that of an unbounded elliptical vortex with uniform 
vorticity. This makes it possible to seek global solutions in the form of plane waves with time- 
periodic wavevector and an amplitude that satisfies a Hill's equation. (Results for this particular 
case have a much broader appeal, however, since an identical Hill's equation arises when the 
stability of more general elliptical vortices is examined locally using the geometric-optics tech- 
nique.—*^) The Floquet problem associated with the Hill's equation is solved asymptotically 
in section 3 in the limit of fast rotation and strong stratification, with the eccentricity of the 
elliptical streamlines assumed of order one. For simplicity, we also make the hydrostatic approx- 
imation assuming that N ^ f and an order-one Burger number. The asymptotic derivation 
is only sketched in section 3, with technical details relegated to Appendix A. The asymptotic 
results are confirmed by direct numerical solution of the Floquet problem in section 4. The 
effect of a finite N/ f is also briefiy examined there. 

2 Formulation 

We consider the stability of a horizontal elliptical flow in a three-dimensional stratified fluid, 
with constant Brunt- Vaisala frequency A^, rotating about the vertical axis at rate //2 > 0. 
The flow's streamfunction, velocity and vorticity are written as 



and note that the flow is anticyclonic for ? = 1, and cyclonic for ^ = — 1. Three dimensionless 
parameters characterise the flow, namely 




U = {ay, —bx, 0) and O = a -|- 6, 




where ab > 0. We define 



? = sgn a = sgn b 





(2.2) 
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which are recognised as the aspect ratio of the eUiptical flow, a Rossby number and the Prandtl 
ratio. We assume that e > 1 without loss of generahty. 

Perturbations to the flow ()2.ip take the form of plane waves with time-dependent wavevector, 
with each fleld written as 

= n(t)e''^W-^, 

where the wavevector k = {k,l,m) satisfies 

k = bl, I = —ak and m = 0, (2-3) 

the overdot denoting differentiation with respect to t. In what follows, we use a dimensionless 
time variable obtained by taking (a6)^^/^ as a reference time. In terms of this variable, the 
solutions to (12. 3p have the simple form 

k = kQCOst, / = — (je/co sin t and m = niQ, (2-4) 



where ko and mo are constant. The stability of ()2.ip depends on the behaviour of the amplitudes 
u{t), v{t), etc. as t — > oo. These satisfy a set of ordinary differential equations with time-periodic 
coefficients. Following McWilliams and Yavneh^, this set can be conveniently reduced to a single 
second-order equation for the amplitude of the vertical component of the vorticity = \{lv — ku). 
Assuming that the perturbation potential vorticity vanishes, this equation reduces to 



•■ 2qklmHe - e^M ■ 1 
C + — -C H 



C = 0, (2.5) 



/fe2 + /2y ^2 p^2 

where = fc^ + /2 _j_ ^2 have omitted the hat on the amplitude C,. Four dimensionless 

parameters appear in this equation: the three fiow parameters ()2.2p . and the initial aspect 
ratio mo/^o of the perturbation. Note that anticyclonic flows (with <;" = !) are susceptible to 
centrifugal instability (e.g. Ref.^^) when the relative vorticity exceeds /, that is, for e(e + e~^) > 
1. Since we focus on the regime e <C 1 we do not this consider this instability further. 

Most of this paper focuses on a limiting case of (j2.5p obtained by making the hydrostatic 
approximation. This assumes that mo ^ kQ and N ^ f while 

^ = # = ^W. (2.6) 

Nko 

This is the regime most relevant to the dynamics of the atmosphere and oceans since the 
condition N ^ f is verified while, as we demonstrate below, the largest growth rates of the 
elliptical instability correspond to ;U = 0(1). The parameter fi can be recognised as the inverse 
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square root of a Burger number; it can be interpreted as the aspect ratio of the perturbation 
scaled by f /N as is natural in rapidly rotating, strongly stratified fluids. 

In the hydrostatic approximation, is approximated by m?, and ()2.5p reduces to 



■■ 2cikl{e-e-^) - 1 



C = 0. (2.7) 



k"^ + P J f^m? 
Using ()2.4p and defining ^ > by 

e2 = l + V^2^ ^2.8) 

we rewrite this equation as 

C - p{t)C + ^ [u^\t) - eq{t) + eV(t)] C = 0. (2.9) 



Here 

^ =1 + U n = l + l + ^^sin^t , 2.10 

can be recognised as the square of the inertia-gravity-wave frequency (non-dimensionalised by 
/). We have also introduced 

V^sin(2t) 1 2e \ , 2(6^ -M) 

K0= ..^ ,2 ■ g(t) = de + e-^+ . and r(t) = ^ . \ (2.11) 

1 + ip"^ sm i \ 1 -|- ^/^^ sm t / 1 + ip"^ sm t 



Equation ()2.9p is a Hill equation, with coefficients that are vr-periodic in t. Its stability is 
determined using the Floquet theory for Hill equations^-. Briefly, if C,{t) = {Ci{t) , C,2{t))'^ is a 
column vector of independent solutions, 

C(i + vr) = MC(t) 

for some constant matrix M. The eigenvalues A of M are then the Floquet multipliers, and two 
fundamental solutions can be found for which 

at) = e'^V(t), (2.12) 

where 

cj = ilogA (2.13) 

vr 

is the Floquet exponent, and (f){t) is vr-periodic. Note that the form of the coefficient of C 
ensures that the two multipliers satisfy A1A2 = 1. Instability occurs when one of these is such 
that |A| > 1 or, equivalently, Reo" > 0. The matrix M is computed by relating C, and C, at two 
times t and t + tt. Here we choose t = — 7r/2 and compute M as 

M= [C(vr/2),C(vr/2)][C(-vr/2),C(-vr/2)]-^ (2.14) 



Our aim is to determine the largest values of the growth rate Reo" for (|2.9p analytically 
in the fast-rotation regime e <C 1, with N S> /, /i = 0(1) and ip = 0{1). In this regime 
()2.9p ressembles the Hill equations with large parameters whose stability has been studied by 
Weinstein and KellerJi using a mapping to parabolic cylinder functions. However, there are 
difficulties in applying their results directly, related to the presence of a first derivative term 
that is singular for the complex values of t such that fc^ + = 0. We have therefore found it 
simpler to develop a different approach, combining WKB analysis with complex-time matching. 
This approach, which has the advantage of demonstrating the link between the instability and 
the Stokes phenomenon— ii^, is described in the next section and in Appendix [A) The analytic 
results obtained there are confirmed and extended to finite N/f in section 4 by solving ()2.9p 
numerically. 



3 WKB analysis 

A WKB solution of the form 



C = yl(t)e^^(*)/^ (3.1) 



can be introduced into ()2.9p . and the (real) functions A{t) and 9{t) can be derived by expansion 
in powers of e. At leading order, we find that 

eo{t) = [ uj{t')dt'. (3.2) 

J-tt/2 

At the next order, we have 

= -I- (3.4) 

We note that Ao(t) is vr-periodic; this is also true for higher-order corrections, so that there 
is no instability to any algebraic order in e: the fundamental solutions are given by C in ()3.ip 
and its complex conjugate, and the Floquet multipliers are simply A = it exp(i^(7r/2)) (taking 
9{—'k/2) = 0). Instability is necessarily an exponentially small effect; furthermore, it can only 
occur for values of the parameters such that A is exponentially close to ±1. This is because the 
condition for instability |A| > 1 requires the two multipliers to be purely real; however, they 
are complex conjugate to all orders in e, a property which persists in the presence of a small 
perturbation in the non-degenerate cases A 7^ ±1. 
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Computations detailed below show that the solution defined by for — 7r/2 < t < —6, 
with <C 1, switches on an exponentially small term as the Stokes line Ret = is crossed. 
Denoting by S the corresponding Stokes multiplier, which is exponentially small in e, this 
implies that the pair of solution 



C = A(t)e'^W/^ and C = ^We"'^^*^/' 
valid for —tt/2 < t < —S (for some e^/^ ^5^1) becomes 



C = Ait) 



and C = Mt) 



(3.5) 



(3.6) 



for 6 < t < 7r/2. Taking = (C, Q"^, we compute the matrix M in (j2.14p . We first have that 
[C(vr/2),C(vr/2)] = 



A(e-^^/^ + 5e'^/^) (-le-^M + i)e-^^/^ + 5(ie-iM + A)e'^/' 



where A, A, 9 and 6 are evaluated at t = 7r/2. Similarly, 

A le-^eA + A 



[C(-vr/2),C(-vr/2)] 



A -ie~^9A + A 



Here A, ^ and 9 can be evaluated at t = 7r/2, as above, since their values at t = ±7r/2 coincide. 
Computing (12.140 gives the simple result 



M 



-le/e 



(3.7) 



Here 9 = 9{'k/2) or, more generally, 9 = 9{'k/2) — 9{—tt/2) which accommodates any convention 
for the arbitrary choice of ^(— 7r/2). 



Suppose now that the parameters are such that 

e^^A = ±(1 + iT) + 0(r2) (3.8) 

for some T G M of a similar order of magnitude as S. The Floquet multipliers obtained from 
()2.14p are then given by 

A = ±(l + V|5p-r2) +0(|5p) and A = ±(l- v'|5|2-r2) +0(|5p). (3.9) 

Clearly, one of these multipliers satsifies |A| > 1, and the flow is unstable, provided that — 15*1 < 
T < \S\, that is, in exponentially narrow instability bands. The corresponding growth rate 
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a 





Figure 1: Parameters a and P governing the maximum growth rates according to p.lOp - p.lip 
as functions of e and /x. 



a = \/\S\'^ — T^/vr + 0(|5p) is maximum at the centre of these bands, for T = 0, and is given 
by 



1^1 



vr 



The computation of S is carried out in Appendix O There we show that 



\S\=e 



-a/e+<;l3 



(3.10) 



(3.11) 



where 



a 



Jo 
Jo 



1 + ^2 — ip'^x'^ 



1 + x2 



e + e ' + 



dx and 



2e 



dx 



(3.12) 
(3.13) 



l-V'^xV ^(l+/i2_^2^2)(l+2;2) 

Here f denotes the Cauchy principal value of the integral, whose integrand is singular at x = 

Figure [T] shows the values of a and /? as functions of e and /i. Some conclusions can be drawn 
from the figure and the examination of the explicit expressions ()3.12p - (l3.13p . First, a — > cxd in 
the limits of both small and large fi; specifically a = 0{fi^^) as fj, ^ and a = 0(log/i) as 
/i — > oo. This suggests, as is confirmed by Figure [U that the largest growth rates are attained 
for fi = 0(1). Thus, the aspect ratio of the perturbations that grow as a result of the elliptical 
instability of vortices should be expected to be the Prandlt ratio: mo/ko = 0{N/f). Second, 
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the behaviour of a for small and large eccentricity is given by 



a ~ ^ (^log^ + l-21og2- -log(l + /i^)J as V'^O, (3.14) 

a ~ ] as w ^ oo. o.lo) 

2^V 

The large- expression ()3.15p can actually be used to estimate a for values of ip as small as 1, 
which makes it very useful. (For /x = 1, for instance, the error in (jS.lSp is 15%, 10% and 5% 
for ij: = 1, 1.5 and 2, respectively.) This expression shows in particular that the largest growth 
rates are attained precisely for /x ~ 1 when is large. Third, the obvious fact that /3 > shows 
that anticyclonic flows (<; = 1) are more unstable than cyclonic flows {q = —1). According to 
(jS.lip . the growth rate in an anticyclonic flow is a factor exp(2/3) larger than the growth rate 
of the corresponding cyclonic flow. Formally, this is an 0(1) factor, but the typical values of (3 
are such that it is numerically very small, so that the instability of cyclones is exceedingly weak 
and probably negligible in most circumstances. Note that because /? is a decreasing function of 
e, the asymmetry between cyclones and anticylones is the largest for small eccentricity. 

The formulas (|3.1Up - (j3.13p give completely explicitly expressions for the maximum growth 
rates of the elliptical instability in terms of the three parameters e, ^ and e (recall that ip = 



Ve^ — 1) . These growth rates are achieved when the three parameters are related in such a way 
that exp(i6'/e) = ±1, that is, 

e = m:e, n = l,2,---. (3.16) 

This condition can be recognised as a resonance condition between the phase of the inertia- 
gravity oscillations and the period of rotation around the elliptical vortex (27r in the dimension- 
less time used here). 

The growth rates can be written more directly in terms of e, /i and e by solving ()3.16p 
peturbatively, with 9 = 9q + eOi + • • •, and and Oi obtained from (|3.2p and (|3.4p . This gives 
the approximate position of the instability bands as well as their width. To leading order, the 
instability bands are centred at values of e and ^ satisfying 



1 r/"^ I 2 /"^ /l _L m2 I ^2^2 

00 = - a/1 + -K^2gin2tdt = - / J ^- J dx = mTe, n = 1, 2, • • • . (3.17) 

^ J-tt/2 jo V 1 — X 

The computation of the correction eOi is more involved. Note that it is in principle needed 
to obtain a leading-order approximation to the growth rate Recr as a function of e and fi. 
This is because the error on a needs to be o(e), which requires to approximate the resonance 
values of e and /z with an o(e) error also. We do not pursue these detailed computations here: 
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since the values of e and fj, satisfying the resonance condition (j3.16p are e-close together, (|3.10p 
provides a useful approximation to the growth rates of the instability without the need to 
locate the resonances accurately. This is demonstrated in the next section where we compare 
the prediction (j3.10p with numerical solutions of the Floquet problem associated with (|2.7p . 

Note that the band width can be deduced directly from the expression (j3.17p for 9q. For 
fixed e and e, for instance, T in (jS.Sp can be written as T = e^^ A/i d/^iO, where is the distance 
between ^ and the resonant values, and the derivative is evaluated at resonance. According to 
()3.9p . the instability-band width is therefore A^ = 2e\S\/d^6 where 9 can be approximated by 



4 Comparison with numerical results 

We have solved the Floquet problem associated with equation (j2.9p for the amplitude C nu- 
merically using Matlab's standard Runge-Kutta solver. The growth rates Re a obtained in this 
manner are compared with the asymptotic estimate (jS.lOp for (Xmax- To emphasise the exponen- 
tial dependence on the inverse Rossby number 1/e, it is convenient to display Re a as a function 
of 1/e for fixed values of fj, and of e. Figure [2] shows the results obtained in case of anticy clonic 
vortices (with ? = 1) for several values of e and fi. Similar results for cyclonic vortices = —1) 
are displayed in Figure [3l 

The figures confirm the validity of our asymptotic estimate. They also suggest that this 
estimate remains useful for moderately small values of e, say e < 1/2. Note that the dimensional 
growth rates are obtained by multiplying a by \/a6 which is related to the relative vorticity 
J7 = a + 6 of the flow by ^/ab = ^}/{e + e^^). As expected from our asymptotics, the growth 
rates in the case of cyclonic flows are exceedingly small for e ^ 1 even for the large eccentricities 
used in Figure O Nonetheless, our results clarify the fact that all elliptical flows are unstable, 
regardless of the sense of the rotation, of its strength, and of the strength of the stratification. 
Note that the match between asymptotic and numerical results for cyclonic flows appears to 
degrade for small e (i.e., large 1/e); this is because the smallness of both the growth rate and 
instability-band width makes the maximum growth rate delicate to estimate numerically. 

The separation between instability bands can be estimated from the asymptotic formula 
(j3.17p : in terms of the varying 1/e used in the flgures, it is given by 
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(a) 



(b) 



Re a 



2.5 3 



(c) 




Re a 



0.1 



2 3 4 



(d) 



0.4- fl 



0.3 



Re a 



0.1 



Figure 2: Growth rates Reo" in anticyclonic flows as functions of the inverse Rossby number 
1/e for (a) e = 1.5, = 1; (b) e = 2, = 1; (c) e = 2, /i = 0.5; and (d) e = 4, ^ = 1. The 
growth rates computed numerically (solid lines) are compared with the asymptotic estimate of 
the maximum growth rates CTmax (dashed lines). 



Evaluating this quantity for the parameters chosen for the figures gives 7 = 0.62, 0.54, 0.31 and 
0.34 for the parameters of Figure [2] (a)--(d), and 7 = 0.18 and 0.12 for the parameters of Figure 
[3] (a) and (b), in good agreement with the numerical results. 

Our derivation of an asymptotic expression for the growth rate makes the hydrostatic ap- 
proximation, which assumes that N ^ f, rriQ S> and fi = fmQ/{NkQ) = 0(1). This assump- 
tion, which could be relaxed, is made because it corresponds to the regime most relevant to 
atmospheric and oceanic flows; it is consistent in the sense that the growth rates obtained are 
maximized for fi = 0(1) and decay rapidly for fi ^ I or fi <^ 1. To test the sensitivity of the 
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(a) 



(b) 



0.08 r 



0.06 



Re a 



Re CTr 



1.5 



1/. 



Figure 3: Growth rates Recr in cyclonic flows as functions of the inverse Rossby number 1/e for 
(a) e = 4, /i = 0.5; and (b) e = 6, = 0.5. The growth rates computed numerically (solid lines) 
are compared with the asymptotic estimate of the inaxiniuni growth rates (^max 

(dashed lines). 



results to the hydrostatic approximation, we have solved the Floquet problem numerically for 
the full, unapproximated equation ()2.5p for moderately large values of A^// and mo/^o- The 
results obtained for /U = 1 and e = 4 are displayed in Figure HI This compares the growth rate 
obtained in the hydrostatic approximation with those obtained for Nj f = rrio/kQ = 3 and 6. 
Except for e > 1, there is relatively little difference between the results: the maximum growth 
rates fall on the same curve, well described by the hydrostatic asymptotics. Of course, the 
location of the instability bands changes depending on S, but this is not significant, since they 
would also change if mo/ko was varied independently of N/ f, as is physically relevant. 

Acknowledgments. The authors acknowledge the support of the EPSRC Network 'Wave- 
Flow Interactions'. JMA is supported by a studentship of the UK Natural Environment Re- 
search Council. 



A Exponential asymptotics 

In this Appendix, we compute the (exponentially small) Stokes multiplier S which quantifies 
the switching on of one branch of the WKB solution by the other (see ()3.5p - p.6p ) through a 
Stokes phenomenon— ii^. The Stokes phenomenon is associated with the presence of complex 
turning points, that is, complex times where lo = 0. From ()2.10p . these are located at 



t„ = isinh — '— — hniT, n = 0, ±1, ±2, ■ 
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Re a 



llll 



1/e 

Figure 4: Effect of the hydrostatic approximation: the growth rate Re a is plotted as a function 
of the inverse Rossby number 1/e for an anticyclonic flow with e = 4 and ^ = 1, in the 
hydrostatic limit N j f oo (solid lines), for N/f = 6 (dashed lines), and for N/f = 3 (dotted 
lines). 

and in- The Stokes phenomenon occurs as t G M crosses one of the Stokes lines joining tn to 
In the interval [— 7r/2, 7r/2] of interest, the only Stokes line is the segment of Re t = joining Iq to 
to- We compute S using matched asymptotics, examining how the solution ( = A(t) exp(i^(t)/e) 
connects to the solution C = A{t)[ex.p{i9{t) /e) + S exp(— i6'(t)/e)] as this segment is crossed.— 
To analyse the behaviour of the solution in the neighbourhood of to, we first note that 



a; ~ ai/V"/^(t - ^o)^/^ where a = y^2(n^Ji?j{TTWTJi^ / fi^ , 



(A.l) 



as t ^ to- We then rescale the evolution equation for ^ near tg by defining the inner variables 

r = e-2/3ai/3(t - to) and Z(r) = ({t). 



(A.2) 



Retaining only the leading order terms, this transforms ()2.9p into the equation 

+ itZ = 0. 



dr2 



Solutions can be written in terms of the Airy functions Ai(e ™/^t) and Bi(e ''^/^r). We claim 
that the solution of interest is 



Z ^ C 



Ai(e-'^/V) + iBi(e-'^/V) 



(A-3) 



We verify that this solution matches A{t) exp(i^(t)/e) to the left of the Stokes line Ret = 0; in 
doing so we find an expression for the constant C. It is convenient to verify the matching along 
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the line phr = — 57r/6; this is an anti-Stokes hne along which the two independent solutions of 
()A.2p have similar orders of magnitudes. Along this line, we can use the asymptotic formulas^ 

Ai(-z) ~ ^_cos(2z3/2/3-^/4), (A.4) 

Bi(-z) _L_sin(2z3/2/3-vr/4), (A.5) 

with z = — exp(— i7r/6)r = exp(5i7r/6)T +00. Thus we have 

On the other hand, using p.3p - (|3.4p the solution = A{t) exp{i9{t)), valid in the outer region 
away from to and to the left of the Stokes line Re t = 0, can be written as 

r ^ ^ /!,,,p(t')df72 ie-\f!,,,[c^(f')-eg(t0/(2c.(f'))]dt' 

/r pit')dt'/2 le-i/r [^(t')-^9{t')/(2c^(t'))] dt' 210-/4^2/3/3 

- (,«)i/6^i/4« - e - e , (A.7) 

after using (jA.ip . Here r_ denotes a contour joining — 7r/2 to to while remaining to the left 
of the Stokes line Ret = 0. Comparing (|A.7p with (|A.6P shows that (" correctly matches Z 
provided that 

^"T^" - e - . (A.8) 

We now match Z with the outer solution valid to the right of the Stokes line Ret = 0. A 
connection formula for Airy functions-^ gives the alternative form 

Z ~ 2Ce'^/3Ai(e-iW6^)^ (A.9) 

for (jA.3p . Carrying out the matching on the anti-Stokes line phr = exp(— ivr/G), we can use 
the asymptotic formula for Ai in (|A.4p to write that 

i7r/24 

Z ^ %— U^^^-I'r^l-I^ + ,e-2--/V^/V3) (A.IO) 

for |r| — > 00. This should be matched with the form ("(t) = A(t)[exp(i0(t)/e) + S'exp(— i0(t)/e)] 
of the solution to the right of the Stokes line. Using (|A.ip gives the asymptotics 

^ _ J_e/-./2P(*')di72 r ie-l/!,/2['^(*')-^9(i')/{2^{t'))]di' ^^g-ie-i/!^/2['^(*')-^'?(*')/(2^(*'))]dt' 

a;i/2 L 



(ea)i/6ri/4 



■^le-i /r+ ['^(*')-f9(i')/{2a;(t'))] dt'^ae-ZVa/VS 



_^^g-i^"'/r+['^{*')-^g(t')/{2a;{t'))]dt'^_2io-/4r3/2/3' 
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where r+ denotes a contour joining — 7r/2 to to- This contour crosses the Stokes hne Ret = 
below the singular point tp of p{t) and q{t), given by tp = isinh"^(l/'(/'). Taking (|X8]) into 
account, the matching with (1A.10|1 gives the two equations 

J^_p(t)dt/2^-if^_q{t)/u{t)dt/2 ^ J^^p{t)dt/2^-iJ^^q{t)/uj{t)dt/2 (All) 
^Jr_ Pi^') dt'/2^,e-^ [^(t')-e<?(t')/(2c^(i'))] dt' ^ ^J^^ p{t') dt'/2^-,e-^ J^^ [a.(t')-eg(t')/(2..(t'))j^i'^2) 

We can now deform the integration contours. The difference /p^ — Jp reduces to an integral 
along a closed contour encircling tp. Computing the corresponding residues using (|2.11|) gives 

Rest p = 1 and Resf — = — i?, 

Taking this into account confirms that (jA.lip is an identity. Using Res tpP = 1 in (|A.12p gives 
the Stokes multiplier as 

g ^ _^^2,e-^f_l^^Ht')~eq{t')/i2u;it'))]dt' ^ 

where the Cauchy principal value integral, denoted by is necessary because q{t) has a pole 
at t = tp. It follows that IS"!, giving the instability growth rate, can be written as 

\S\ = e-"/^+^^, 

where the two constants 



a = -2i [° uj{t) dt and /? = -i/ ° '^^^^^ 
Jo Jo 



dt 



a;(t) 

are real, positive and independent of <;^. Using (I2.10p and (12. lip , they can be given the more 
explicit forms 



2 



a = — \/ 1 — sinh u+ u? dn, 

Jo 



and 



.sinh \./T^/^) / 2e \ dn 

I3 = ii<.j (e + e ^+ * 



1 - V^sinh^ny a/1 - sinh^ n + /i^ ' 
and further transformed into the convenient expressions ()3.12p - (l3.13p . 
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